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roots of In + 1 , the algebraic sum of which roots is one, and dimi- 
nishing the middle difference by theorem Q until it reaches a number 
nearest to half the sum of the first and third differences. The 
difference between 2n -\- 1 and the number thus obtained will be the 
double of a triangular number = 2T. By the next step, the extreme 
differences are reduced until they are of the form m, m-\-\ ; and the 
difference between 2ra + l — 2T and the number thus obtained will 
again be the double of a triangular number = 2T>. The differences 
last obtained give the double of a triangular number + 1 — 2T" + 1 . 
So that we find 2b + 1 = 2T + 2T' + 2T" + 1 . Consequently »= the 
sum of three triangular numbers, if all the three operations be 
necessary; if not, to two or one triangular number only. 



II. The first part of a paper " On a Class of Differential Equa- 
tions, including those which occur in Dynamical Problems." 
By W. F. Donkin, M.A., F.R.S., F.R.A.S., Savilian Pro- 
fessor of Astronomy in the University of Oxford. 

This paper is intended to contain a discussion of some properties 
of a class of simultaneous differential equations of the first order, 
including as a particular case the form (which again includes the 
dynamical equations), 

dt/i dXi 

where x x ... x u> y l ...y n are two sets of n variables each, and accents 
denote total differentiation with respect to the independent variable t ; 
Z being any function of x, &c, y v &c, which may also contain t ex- 
plicitly. The part now laid before the Society is limited to the 
consideration of the above form. 

After deducing from known properties of functional determinants 
a general theorem to be used afterwards, the author establishes the 
following propositions. 

If x 1 ... x n be n variables connected with n other variables y 1 ...y n 

by n equations of the form y ; = — - (X being a given function of 
x,...x„); then the equations obtained by solving these algebrai- 



eally, so as to express cc l ... x„ in terms of y , . , . y,„ will also be of the 

dY 
form Xi= — ; where Y is a function of y, ... y„, which may be de- 

fined by the equation 

Y=-(X) + Oi)y>+ ... +(Oy„, 

in which the brackets indicate that the terms within them are to be 
expressed as functions of y a ...y n . Moreover, if p be any other 
quantity contained explicitly in X (besides the variables x l ... x„), 
the following relation will subsist ; namely, 

dX dY_. Q 
dp dp 

the differentiation in each case being performed only so far asp 
appears explicitly in the function. 

It is then shown that if X contain explicitly, besides x x ... x„, the 
n constants a,, « 2 , ... «„, and the variable t, and if the 2« variables 
a?j ... x n , y, ... y,„ be determined as functions of t by the system of 
2m equations, 

ttX flX t /TT \ 

dxi " dcti 

where b x ... b n are n other constants, the elimination of the 2ra con- 
stants from these equations and their differentials with respect to t, 
leads to the system of differential equations (I.), if for Z be put the 

result of substituting: in — ■ — the values of the In constants in 
5 dt 

terms of the variables. The equations expressing the 2« constants 

in terms of the variables may be considered as the 2» integrals of 

the system (I.). 

The author employs the symbol [jo, q] in a sense similar to that 

in which Poisson and others have employed (p, q), namely, as an 

abbreviation for 2;( -£- J- -- — ) ; and he shows that if p, q 

x rfy ; dx t dx t dy/ 

represent any two of the 2» constants «, &c, J, &c, then \_p, j] is 

either =1 or =0, according as^i, q are a conjugate pair a t , b it or not. 

Next it is shown that if «,, « 2 , ...«„ represent any functions of 2nn 

variables x s ... x n , y, ...y n , satisfying identically the - *■ ~ ' con- 
ditions [a,, aj]=0, then if by means of the n given equations ex- 



pressing a lt &c. in terms of the variables, the set y v y 2 ... y„ be 
expressed as functions of x l ... x n , a, ...«„, the - ' relations 

~ = -^-i will be identically satisfied ; in other words, the expres- 

sion for y, ... y n will be the partial differential coefficients of a func- 
tion of*! ... x n . 

Hence it easily follows, that if any n integrals a t ...«„ofthe 
system (I.) be given, which satisfy the conditions [« f , «J = 0, a 
" Principal Function " X can always be found, from which the re- 
maining integrals of the system may be deduced by means of the 
second set of equations (II.). 

The relation in which these investigations stand to the discovery 
of Sir W. R. Hamilton (as improved and completed by Jacobi) is 
pointed out. And it is shown that the system of n differential equa- 
tions of the second order 

ZdW\'_dW 
^dx'i' dx t 

(to which Lagrange had reduced the dynamical equations, and which 
Sir W. Hamilton had transformed into the system (I.) by a process 
depending upon the circumstance that, in dynamical problems, W 
contains x\, x\, ... x' % only in the form of a homogeneous function) 
may, by means of the theorems established at the beginning of the 
paper, be reduced to the form (I.) without assuming anything as to 
the form of W, which may be any function whatever of x x ... #„, 
w\ ... x' n , and f. 

The 2ra integrals of the system (I.), obtained in the way above 
explained, being shown to satisfy the conditions 

[«,, JJ = 1 , [«j, a,] = [a*, JJ = [*„ 5,] = 0, 
it is proposed to call them " normal integrals," and the constants 
«! &c, b t &c. " normal elements," any pair a it b t being called con- 
Jugate. 

In the second section, the author gives a simplified demonstration 
of Poisson's theorem (extended to the general system (I.)), that if 
f, g be any two integrals, [/, g~\ is constant. The preceding prin- 
ciples are then exemplified by application to the problems of the 
motion of a material point under the action of a central force, and 
the rotation of a solid body about a fixed point. 



In each case three integrals, c v c 2 , c 3 , are taken, satisfying the 
three conditions [c 2 , c 3 ]=0, [e 3 , c,] = 0, [e„ c^^O; the first being 
the integral of vis viva, and the other two being derived from the 
integrals expressing the conservation of areas. In the former pro- 
blem the " principal function " is then found with great ease, and the 
remaining integrals deduced. The set of " normal elements " thus 
obtained coincide with those given by Jacobi (in a memoir in Crelle's 
Journal, vol.xvii.). In the problem of rotation, the algebraical solution 
of the three assumed integrals for y x , y 2 , y 3 depends upon that of an 
equation of the fourth degree. It is therefore impracticable to 
exhibit the principal function in an explicit form. In this respect 
the result arrived at resembles that obtained by Mr. Cayley in a 
totally different way ; Mr. Cayley having shown that the solution of 
the problem is reducible to quadratures, assuming the algebraical 
solution of a certain system of equations of the same form as those 
to which the author of the present investigation is conducted. (Camb. 
and Dub. Math. Journ. vol. i. p. 172.) 

Methods are then indicated by which, when one system of "normal 
elements " is given, other systems may be found. 

The practical value of " normal solutions " of the system (I.) de- 
pends chiefly upon the simplicity of the corresponding formulae for 
the variation of elements, the theory of which is intended to form 
part of the subject of the following sections. 



